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Abstract 

We study the gauge-invariant gaussian ansatz for the vacuum wave functional and 
show that it potentially posseses many desirable features of the Yang-Mills theory, like 
asymptotic freedom, mass generation through the transmutation of dimensions and a 
linear potential between static quarks. We point out that these (and other) features 
can be studied in a systematic way by combining perturbative and 1/n expansions. 
Contrary to the euclidean approach, confinement can be easily formulated and easily 
built in, if not derived, in the variational Schrodinger approach. 
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1 Introduction 

We shall consider the pure Yang-Mills theory, else called gluodynamics, with no dynamical 
quarks. It is widely believed that in this theory the potential between probe static quarks 
is infinitely rising linearly with the quark separation. If proved correct, it will be one of the 
most remarkable findings in physics, and one would need to explain it from different angles. 
The present paper is an attempt to understand the physics of confinement in the framework 
of the Schrodinger picture. 

The use of the Schrodinger equation to discuss the ground state of gluodynamics as well 
as the two source quarks on top of the ground state, is not new. To mention a few early 
works, see refs.[|T], 0, [| |J. More recently, the interest to the Schrodinger equation has been 
revived by Kogan and Kovner who suggested a variational vacuum wave functional in a 
gaussian form averaged over all gauge transformations . We shall use this form of the trial 
wave functional making a slight but essential generalization of the ansatz suggested in ||. 
This work has been influenced also by a recent paper by Zarembo || though we differ from 
that reference in several aspects. 



2 Schrodinger equation 

The spatial components of the Yang-Mills potential A\ (x) play the role of coordinates, the 
electric field strength Ef(x) = Af(x) plays the role of momenta. The classical energy of the 
Yang-Mills field is 



H dass = J d 3 x 

where 

Bt = tiAdjAl + {l/2)r hc A)AD (2) 

is the magnetic field strength. Upon quantization the electric field is replaced by variational 
derivative, Ef(x) — > —ig 2 5/5A^(x), if one uses the 'coordinate representation' for the wave 
functional. The functional Schrodinger equation for the wave functional takes the 

form 

m[M = /^{-y^^ + ^(^)) 2 }^[^ =£*IM ( 3 ) 

where S is the eigenenergy of the state in question. 

The Schrodinger equation (|3]) should be supplemented by the gauge condition. For a 
state without external colour sources the condition is that the wave functional ^f[Ai] should 
be gauge-invariant, 

*[Af]=*[Al (4) 

where 

Af = S j AS + tS j d,S (5) 
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is the gauge transformation. For a state with an external static quark at z\ and an antiquark 
at z 2 the wave functional ^[Al z ii ^2] should satisfy the gauge condition 



^[Af; Zl , z 2 \ = [&(zi)] a VfiAn z u z 2 \ [S(z 2 % . 



(6) 



If the probe sources belong to another representation the gauge-tranforming matrices S(z\ t2 ) 
should be taken in appropriate representation. 



3 Ultraviolet regularization 

Strictly speaking, the hamiltonian (||) is senseless since any calculations with it are plagued 
by ultraviolet divergences. In order to give meaning to the Schrodinger equation (0) one has 
to regularize both the kinetic and potential energy parts of the hamiltonian. Examples of 
ultraviolet regularization are provided by the lattice hamiltonian [0 or by the 'heat kernel' 
method ||. 

The requirements on possible ultraviolet regularizations are rather stringent. In addition 
to the apparent requirement that the regularization should be gauge invariant, there is 
a less trivial requirement that the regularized hamiltonian, together with the momentum 
and angular momentum operators, should satisfy the Poincare algebra. Otherwise, there 
would be various inconsistencies, like that the energy of a glueball state would not obey the 
relativistic formula E 2 = m 2 + p 2 . To my knowledge, a consistent regularization scheme has 
not been yet introduced. 

Our aim is, however, the vacuum state with and without static sources, therefore any of 
the two mentioned regularization schemes would suit us here. In fact, one has to solve the 
regularized Schrodinger equation 



H re9 (g, M)*{g, M) = £{g, M)*{g, M) (7) 

where g is the bare coupling constant given at certain large normalization scale M which 
one should eventually put to infinity. [For example, M can be viewed as the inverse lattice 
spacing in case of the lattice regularization]. 

Let us consider the ground state or vacuum of the theory. Its energy £ is by far domi- 
nated by the energy of the normal zero-point oscillations of the gauge field, which diverges 
quartically and, being regularized, is proportional to M 4 . From the perturbation theory one 
infers 

£^(g, M) = 2(N 2 - 1) V f-1 + pertuAative 
v ' y 6 9 corrections 



2(iV c 2 -1)VJ- 



d 3 k \k\ perturbative 



(27r) 3 2 corrections 

= V M 4 (coo + c 01 g 2 + c 02 g 4 + ...). (8) 

It is not the full truth, however. We expect also non-perturbative contributions to the vacuum 
energy determined by the dimensional parameter A obtained through the transmutation of 
dimensions, 



3 



where b is the one- loop and b' is the two- loop Gell-Mann-Low coefficient. Contributions 
proportional to powers of A cannot be obtained in any order of the perturbation theory. 
However, there are no reasons why such contributions should be absent. Therefore, we write 
for the vacuum energy a general double expansion: 

g( ^ M) = M 4 (c 00 + c 01 g 2 + c 02 s 4 + ...)+ M 3 A(c 10 + c ll9 2 + ...) + M 2 A 2 c 20 + ... + MA 3 c 30 + ... 

+ A 4 (c 40 + C4 i3 2 + •••)• (10) 

Meanwhile, all physical observables, like the string tension in the static potential, glueball 
masses, gluon condensate, topological susceptibility, etc., should arise as renorm- invariant 
combinations of the bare coupling g and the ultraviolet cutoff M, i.e. should be expressed 
trough A only, in appropriate powers. This makes a straightforward use of the variational 
principle to extract physical observables almost a hopeless task: before getting to physically 
interesting quantities encoded in the last term in eq. (pTOj) one has to minimize the first four 
terms which are by far the larger! 

There is an interesting possibility to circumvent this difficulty, at least for the vacuum 
state. Let us assume that we have succeeded in building the regularized operators for the full 
stress-energy tensor, 6^ v = [F^ a F£ Q — {\/A)g tlv F^ a F^ \/g 2 , in accordance with the Poincare 
algebra, the hamiltonian being 7i = JcPxOqq = Jd 3 x (E 2 + B 2 )/2g 2 . The vacuum energy is 
£ = V(9 00 ) where (...) denotes averaging over the vacuum wave functional. 

Because of the Lorentz invariance of the vacuum one expects (6*00) = —(#11) = —(#22) = 
— (#33) and hence (6*00) = (1/4) (0^). The last quantity is naively zero (following from the 
above definition of 9^ u ), however under proper regularization it should be equal, through 
trace anomaly, to 



B W)(FP 2 _ h B 2 -E 
4^ 4 32vr 2 16n" ' 



[We have neglected higher- loop contributions to the Gell-Mann-Low function f3(g 2 ) which 
can always be done if the ultraviolet cutoff M is chosen large enough and hence the bare 
coupling g small enough.] The quantity {(F^) 2 ) is known as the gluon condensate appearing 
in the ITEP sum rules and has to be renorm- invariant and hence proportional to A 4 . 

Thus if one wants to extract directly the non-perturbative physics from the variational 
approach one should rather maximize the gluon condensate (B 2 — E 2 ) than minimize the 
energy density (B 2 + E 2 )/2g 2 |. 

In the latter quantity the kinetic and potential energies of the zero-point gluon oscilla- 
tions add up and diverge quartically, see eq. (§D; in the former quantity they cancel each 
other. In the zero order of the perturbation theory this cancellation is exact due to the 

2 This is what one actualy does in building the euclidean instanton vacuum from the variational principle 
||, where one is also interested in finding the shift of the vacuum energy in respect to the perturbative one. 
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virial theorem for harmonic oscillators. Under a proper regularization of the operators, re- 
specting the 50(4, 2) generalization of the Poincare algebra (including the dilatation and 
conformal operators), this cancellation should hold to any order of the perturbation the- 
ory, so that the properly defined gluon condensate gets a non-zero contribution only from 
the non-perturbative physics which is characterized by being non-analytic in the coupling 
constant, see eq. (|9|). 

Unfortunately, the prescription to maximize (B 2 — E 2 ) is a priori not well-defined either 
as one can always invent a crazy trial functional making the quantity as large as one likes. 

Despite these pessimistic remarks concerning the foundations for the use of the variational 
principle, it may still serve as a useful guide to understand non-perturbative physics in a 
language different from the usual euclidean approach. Being unable to formulate the problem 
in a mathematically unequivocal way, we shall henceforth proceed in a somewhat sloppy 
manner: we shall not regularize the hamiltonian but simply cut the divergent momenta 
integrals 'by hands' at certain large momentum M ^> A. We shall not go beyond the 
one- loop calculations, so this barbarous procedure will be sufficient for our needs. In fact, 
we shall be able to separate the perturbative and nonperturbative physics using a specific 
variational ansatz, and minimizing the nonperturbative vacuum energy will make sense, see 
section 8. 



4 Kogan— Kovner variational ansatz and its generaliza- 
tion 

Kogan and Kovner || have suggested a gaussian trial wave functional for the vacuum state, 
averaged over all gauge transformations of the 'coordinates' Ai : 

9[A]=J DU(x) exp L±Q[A V ]\ (12) 

where Af = WA4U + iWdiU is the gauge-transformed Yang-Mills potential and Q[A] is a 
quadratic functional (hence the notation): 

Q[A]= J jd 3 xd 3 yK ij (x-y)TrA i ^)My)- ( 13 ) 

The kernel Kij(x — y) is the only variational (or trial) function of the ansatz. Its Fourier 
transform has the general structure compatible with space isotropy: 

K l ,(p) = a(\p\)-5 lJ +b(\p\)-^. (14) 

The two scalar functions a(p) and b(p) are the variational 'parameters' of the ansatz. As 
we shall see below, minimization of the vacuum energy in the leading order of the perturba- 
tion theory (the coefficient Cqq in the notations of eq. fllOD ) leads to the following result: 

a pert {p) = \p\, b pert {p) = arbitrary. (15) 

In ref. |5[] the longitudinal part of the kernel b(p) has been put to zero (without compelling 
reasons). In ref. |J it has been shown that, because of the averaging over gauges in the wave 
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functional (0), certain forms of "actions" in the exponent of *&[A] are physically equivalent; 
in particular, the longitudinal part of the quadratic form Q[A] can be removed by choosing 
another "action" in the exponent, containing higher powers of A4. Having no objections to 
this statement in general we note, however, that if one resctricts oneself to purely quadratic 



"actions" of the form (|H|) the longitudinal part of Q[Aj\ becomes physically meaningful. 
Moreover, we shall show below that by suitably choosing the longitudinal function b(p) 
in (nf) one can reproduce the correct (H/3)iV c coefficient in the one- loop Gell-Mann-Low 



function. In ref. the value (12/3)iV c has been obtained instead for that coefficient, with 
b(p) set to zero. 

According to the variational principle the vacuum energy satisfies the inequality 

£ > (16) 

One has to minimize the r.h.s. in the free parameters / functions of the chosen ansatz, in 
our case, in the two scalar functions a(p) and b(p). 

5 Evaluation of the norm M = (^l 1 ^) 

The norm M standing in the denominator of eq. fllCf ) is understood as the functional integral, 

J\f = J DA m*[A]m[A] 

= jDAjDU 1 J DU 2 exp {-^Q[A^\ - jQ[A^]\ • (17) 

One first integrates over the gauge potential A by saddle-point method; since Qt^f 1 ' 2 ] 
are quadratic in A the saddle-point integration is exact. 
We have 

^Lj = j d 3 y K^x - y) 2 Tr tf(x)t a U(x) [u^y)A j (y)U(y) + iU\y)d j U{y)] (18) 

where t a are SU(N C ) generators in the fundamental representation normalized to Tr t a t b = 
(l/2)5 ob . The saddle-point equation for A4 is 

6Qm 6Q[A^] = 

5A?{x) 5A?{x) ' 1 ; 

The solution of this equation gives the saddle-point value of Af. 

Ai= l - (u x diUl + U 2 diUi) + corrections. (20) 

In the one-loop approximation it is sufficient to neglect corrections to this value as they 
prove to be actually of the order of g 2 , see below. This saddle-point value becomes exact if 
the kernel K(x — y) is a delta function considered as an example in ref. B. In ref. H more 
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complicated technique has been involved to solve the saddle-point equation than seems to 
be necessary for our goals here. 

It is amusing that 'on the diagonal' where the gauge transformations in \1/ and the con- 
jugate \I>* coincide (U\ = U 2 ) the saddle-point value A\ (that is where |\I/L4]| 2 is mostly 
concentrated) is a pure gauge potential. That is consistent with the euclidean instanton- 
vacuum picture where the system spends much of the 'time' in a pure gauge state. 

Actually, the introduction of two gauge transformations U\ and U2 is abundant as one 
can always absorb one of them into a shift of A{ since the integration measure DAi is gauge 
invariant. Renaming the integration variable A^ 1 — > Ai so that A^ 2 becomes Af where U is 



the 'relative' gauge transformation, U = U\U%, one can rewrite the norm as 



M = JDU jDAexp 
and the saddle-point field as 



1 



1 



Q[A] - - 2 Q\A 



(21) 



A4 = —Udilft + corrections. 
2 



(22) 



Generally speaking, this is not a pure gauge potential; it becomes a pure gauge only when 
the relative gauge transformation U does not differ much from unity. 

Let us introduce the shifted field A' counted from the saddle-point value (|22|) , i.e. let us 
put Ai = (i/2)UdiW + A\ and expand Q[A] + Q[A U ] in eq. (0) in powers of A'. 

The term quadratic in A' is 



d 3 xd 3 y Kij(x - y) {2Tr A'^A'^y) + Tr \xj{y)U\x) 



A&x) \U(x)U\y)-l] A'(y) 



+Tr [U(y)U\x) - l] A'^A'^y) + TrA^x) [u{x)U\y) - l] A'^y)} . (23) 
The term linear in A' is 



d 3 xd 3 yK lJ {x-y)(-t)TiA' i {x){[u{x)U\y) - l] \u{y)d j U\y)] [u(y)U\x 



+ 



U{x)U\y) 



U{y)d j U\y)\ + [U(y)d 3 U\y)\ [U{y)U\x 
Finally, the term not containg at all A' is 

J j d 3 xd 3 y Kijix - y) {^Tr d i U\x)d j U{y) 



(24) 



+ h:id i U\x)d J U(y) [u\y)U{x) 



-Tr d i U{y)d j U\x) [U{x)U\y) 



(25) 



We see that had Kij(x—y) been a (^-function one could replace U(x)W(y) by unity matrix, 
so that the linear term in A' ( |24] ) would vanish. In the same approximation one would be 
left only with the first terms in eqs.(|23|, |25| ). In fact this approximation to Q[A] + Q[A C/ ] is 
correct for any function Kij(x — y) if we restrict ourselves to the one-loop approximation, 
i.e. neglect 0(g 2 ) terms in calculating the norm M . 

To show that we first notice that if we, indeed, neglect the linear term ( pip and leave 
only the first term in eq. (|23|) and the first term in eq. (p5|), the propagator of the A[ field 
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is of the order of g 2 , and the propagator of the U field is also of the order of g 2 . The last 
statement needs to be made more precise as we, generally speaking, deal with a non-linear 
field: that will be done below. What can immediately be seen is that the fluctuations of 
the gradients of the U field are of the order of g 2 , hence the characteristic square bracket 
[U(x)W(y) — 1] is of that order, too, since it is non-zero only when one expands it in the 
gradients of the U field. 

To check that we have correctly neglected all the terms we expand exp(— Q[A]/g 2 — 
Q[A u ]/g 2 ) in powers of the presumably small terms. For example, let us consider the third 
and fourth terms of eq. ( P3| ) expanded to the linear order. There is a factor 1/g 2 coming 
with the expansion, a factor g 2 arising from the contraction of the two A' fields and an extra 
factor g 2 owing to the fluctuations of the square brackets, as explained above. We obtain 
thus a g 2 correction, as compared to the main contribution to the norm which is 0(1). 

Let us consider the expansion in the linear term (p4|) . Expanding to the first order we 
get an identical zero (because it is odd in A'), so we have to consider the second term in 
the expansion of the exponent. We have (1/g 2 ) 2 from the expansion to the second order, g 2 
from the A' propagator and at least g 4 from averaging over the C/'s, that is again a 0(g 2 ) 
correction. The same is true for the expansion of the last two terms in eq. fl25|). 

We conclude that if we are willing to neglect terms of the order of g 2 (i.e. restricting 
ourselves to the one- loop approximation) we can use the saddle-point field given by eqs.(|20|, 
^) and approximate the norm by a factorized expression: 

M ps Ma • Mu, (26) 



M A =jDA'exp^~ IJd 3 xd 3 yK ij (x-y)TrA' i (x)A' j (y)^ 



(27) 



Mu= jDUexpi--^ JJd 3 xd 3 yK ij (x-y)Tid i U\x)d j U(y)\. (28) 

We see that the norm is factorized into independent integrals over A' and over U. The 
former is a trivial gaussian integral while the latter is non-trivial since the field U is subject 
to the unitarity constraint. 

6 Evaluation of the energy numerator (^|7i|^) 

We now proceed to evaluating the numerator of the vacuum energy fll6"l) . We have 

<*,*,*> = joujDAj^ {IH^fg + WW)'} exp {-1(QW + QW)} 

(29) 

The exponential factors are the same as for the norm integral, so we have to analyze the 
prefactor which is new. Using eq. (18) we get 

/ d3x 5A^(I)l Q A [ fx) = III d " xd "v d " zK ^ x - V) K ^ X ~ z ) 2 Tr [u\x)A 3 (y)U(x)A u k (z)_ . 

(30) 
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We next expand this expression in the shifted field A' around the approximate saddle 
point substituting A$ = (i/2)UdiU' + A' i . To save space, we write down the expansion of the 
last factor in eq. ( pCf ) only, the trace. 

The term quadratic in A' is 



TxA'Jy)A' k (z)+Tr \u(z)U\ 



x 



A'Ay) \u{x)U\z) 



+ Tr 



U(z)U\x) - l] A'^A^z) + Tr Aj(y) \u{x)U\z) - l] A' k (z) 



(31) 



The term linear in A' is 



^Tr A' k (z) {U(z)W(x)U(y)d j U\y)U(x)W(z) - U(x)W(y)U(y)d j W(y)U(y)lP(x)} . (32) 



At y — x — z this term is zero implying that in the general case it can be identically rewritten 
a la eq. (^4[) using expressions of the type U(z)W(x) — 1 . However, it leads to a rather 
lengthy formula, and we do not present it here. 
Finally, the term without the A' field is 

-i [lid j U ] {y)d k U{z) + Tr \u\x)U{y) - l] dp\y) \u{x)U\z) - l] d k U(z) 
+Tr \u\x)U{y) - l] djtf (y)d k U(z) + Tr djU^y) [u(x)U^z) - l] d k U(z)} . (33) 



We have now to integrate the sum of eqs.([H], [32], |33|) over the fields A' and U with the 
exponential weight given by eqs.(|27], ^8|). The leading 0(1) contributions come from the 
first term in eq. (^T|) and from the first term in eq. (|3"B|). Eq. ( |3"2"D can contribute only in 
combination with the linear term from the exponent ([24]) expanded to the first order: it 
gives a g 2 correction according to the power counting rules of the previous section. By the 
same counting rules one can neglect all other terms in eqs. (|3~T|, |33|) . 

A similar scenario takes place when we insert A4 = (i/2)UdiW + A[ into the magnetic 
energy B 2 /2g 2 . The leading 0(1) contribution arises only from the 'abelian' piece made of 
the A' field, namely from Tr {diA'- — djA'i) 2 /2g 2 . All other terms, including the Skyrme-type 
term made solely of the U field, give 0(g 2 ) corrections to the vacuum energy. Disregarding 
systematically all such terms we finally obtain: 

« jDujDAexp^-^ j Jd 3 xd 3 y K^x-y) 



2Tr Mx)Mv) + ^TTd l W(x)d j U(y) 
x ^ {/// d 3 xd 3 yd 3 zK tJ (x - y)K ik (x - z) [Tr A j {y)A k {z) - ^Tr djU^y) d k U(z) 



+±Jd 3 xTr(d i A j -d j A i f} 
(we have deleted the primes in the notation of the A field). 



(34) 
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To get the variational estimate for the vacuum energy we need now to divide (|34]) by the 
norm (26). Since the norm has a factorized form the vacuum energy is a sum of two terms: 



£ > 



(V\H\V) 



contribution 
from A field 



+ 



contribution 
from U field 



(35) 



where 



contribution 
from A field 



JDAexp [...] 



9 fDAexp [-f 2 ffffix d 3 y K tJ (x - y) TrMtfMv) 
d 3 xd 3 yd 3 zK i:i (x - y)K lk {x - z)Ti A 3 (y)A k (z) + \jd 3 x Tr (<9 4 A, - d^fX , (36) 



contribution \ 1 JDU exp [...] JJJdx dy dzKij(x — y)K ik (x — z)Ti djU^ (y)d k U(z) 



\ from U field J Ag 2 JDU exp [-^ Jfdx dy K^x - y)Tr d i W(x)d j U(yj 

(37) 

Here exp[...] in both expressions mean the same as in the corresponding denominators. The 
first term in eq. ( |36"D is the contribution of the A, field to the electric energy while the second 
is the contribution to the magnetic energy. Notice that in the one-loop approximation the 
U field contributes only to the electric energy. 

The gaussian integrals over A, can be easily done in the momentum space using the 
propagator 

(A$(p)A*(-p)) = £fi*K?(p). (38) 
Inversing the general expression (|14]) for K^ip) we have 



1 L PiPj Hp) 

a(p) \ 13 p 2 a(p) + b{p) t 
We are now in a position to express the contribution of the Ai field to the electric and 



VW-^K-^tt^tI- (»> 



magnetic energy of the vacuum through the trial functions a(p) and b(p). Using eqs.(^, [38|, 
j3^ ) we get 

( magn6tiC )=i^l = N '- 1 l^Lm- . - p.p.) 9 -KrM p ) 
{ energy density J 2g 2 2g 2 J (2tt) 3 {P %1 PiPj) 2 13 W 

N 2 -lfd 3 p p 2 



(27r) 3 a(p)' 

electric \ (E 2 ) N 2 - 1 r d?p T ^ T ^ .# 2 , 
energy density J = ^ = Jfr)* ^ * ^ + ® 

N 2 - 1 r d 3 p 3a(p) + 6(p) 



(27T) 



+ (37). (41) 
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To write down the contribution of the U field to the (electric) energy density in a more 
compact way let us introduce new kernels, 

L(x-y) = d xi d yj K ij (x-y), L(p) = PiPjK^p) = p 2 [a(p) + b(p)] (42) 

and 



S(y -z)= f d?xd yj K l3 {x - y) d zk K ik (x - z), S(p) = p 2 [a(p) + b(p)] 2 = (43) 
J p 



Eq. fl57|) can be compactly rewritten through the variational derivative of the norm, 

contribution \ 1 f d 3 p 5 hgMu (2ir) 3 . 

l S (P) F773 TJ-i ( 44 ) 



where 



from U field I 2J (2vr) 3 5L{p) V ' 



Nu = jDUexp -^JJ?PL(p)Ttrf(p)U(-p) 



(45) 



2gV (2tt) 

In the next two sections we shall study the contribution of the U field to the electric 
energy density in two regimes: perturbative, which will serve as a pedagogical example, and 
non-perturbative. 

7 Perturbative regime 

In the perturbative regime one assumes that the integration over the unitary field U in 
eqs. (|28| , |45|) is concentrated in the vicinity of unity matrices. Hence we can expand U as 

U = exp(i(f) a t a ) w 1 + i(j) a t a , U^l- i(j) a t a , DU = D(j) a . (46) 

In this approximation the action for the N 2 — 1 fields <p a becomes quadratic, leading to 
the propagator 

?n 2 t) ab ?n 2 fi ab 

mm . p)) = A_ = ?R A_ . (47) 
The addition to the electric energy density arising from the U field is thus 

perturbative contribution \ N 2 — If d?p 2g 2 

:S(p) 



from the U field J 8g 2 J (2tt) 3 L(p) 

N 2 - 1 r d 3 p L(p) N 2 - 1 r d 3 p a(p) + b(p) 



2 J (2tt) 3 2p 2 2 J (2tt) 3 

Naturally, this result comes also from the general formula (|4| 
Adding this result to eq. fl4T| ) we obtain 



(48) 
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electric 
energy density 



(E 2 ) N 2 — If d 3 p 



Mp)- 



(49) 



2g 2 2 J (2tt) 3 

Notice the cancellation of the longitudinal piece associated with the trial function b(p). This 
is as it should be in the perturbative regime. 

We now add up the electric (fi9|) and the magnetic fl40|) parts of the vacuum energy and 
get the variational estimate for the vacuum energy: | 



V ~ 2 



1 f d?p 



a(p) + 



a(p) 



(50) 



This expression should be varied in respect to the trial function a(p). The r.h.s. of eq. (|50| 
has the minimum at 



a(p) = \p\, b(p) = arbitrary, 
the vacuum energy at the minimum being 



d 3 p \p\ 



(51) 



(52) 



(2tt) 3 2 ' 

This is, of course, the correct energy of the zero-point oscillations of N 2 — 1 free gluons with 
two physical polarizations, each mode with momentum p carrying the energy \p\. Note that 
at the minimum the virial theorem for harmonic oscillators is satisfied, since (E 2 ) = (B 2 ). 

These nice results do not come unexpected since the gaussian wave functional used here 
is exact for the free fields. We have only shown that its exact form given by eq. (|5"ID can be 
found from the variational principle. 



8 Key point: non-zero v.e.v. of the Lagrange multi- 
plier 

~~]) can be parametrized 



Let us consider the SU(2) gauge group. The unitary matrix U in eq. 
by quaternions, 



U 



a;, 



±ir,l , 



ii' 



(53) 



The last condition imposed on four real fields u a ensures that U is a matrix from SU(2). 
The functional integral over SU (2) matrices U can be written as 



JDU = J Du a \{5 (u 2 a (x) - l) = jDu a jD\exp —— J d 3 x(Xu 2 a - A) 



(54) 



where we have introduced an auxiliary integration over 'Lagrange multiplier' field A (a;) to 
ensure that the unitarity condition is satisfied at all space points. The factor 1/g 2 is inserted 
in eq. (H) for future convenience. The integration measure (53) is invariant under left and 
right shifts along the group, therefore it belongs to the class of the Haar measures. 



3 Actually, eq. © can be found in rcf . |s| . 
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The dependence of the energy numerator ^) and of the denominator (^l^) on the 

U field enters through the expression (see eq. (Wft)) 



Therefore, the norm 



Trd i U^(x)d j U(y) = 2d i u a (x)d j u a (y). 
or fl45D reads 



(55) 



N v = jDXexp( — Jd 3 x XjjDu a exp j— 2 



d xd yL(x — y)u a (x)u Q (y) + / d xXu 



(56) 

where we have used the new kernel L as defined by eq. (|2). The inner integral is a gaussian 
integral over four real fields u a with the Lagrange multiplier A playing the role of the mass 
term. 

We have learned from the previous section that a(p) = \p\ at large momenta. The 
behaviour of the trial function b(p) is not defined by the leading-order perturbation theory, 
however one can reasonably assume that it can be also ~ \p\. Therefore, the asymptotics of 
the 'kinetic energy' term in eq. (KEf) is 



L(p)^c\p\\ c=l + b, (57) 

where b is the (unknown) numerical coefficient in b(p) — > b\p\ at large momenta. 

Although the number of fields is only four (n = 4) we shall treat it as a formal parameter. 
The gaussian integral over n real fields u a in eq. fl56|) can be formally written as 



J Du a exp {...} = exp 



n 



-Splog(L + A) 



{hi 



where Sp stands for the functional trace. Let us divide the Lagrange multiplier field \{x) 
into a point-independent part /i 3 and a space- variable part X'(x): 



X( x ) =/j, 3 + X'(x), ld 3 xX'(x) = 0. 
The functional trace in fl58T) can be rewritten as 



(59) 



Splog(L + A) = V 



d 3 p 
(2^)3 



log L(p) + /i s 



' ' f/ ' g -A'(g)n(g)A'(-g) + 0(A' 3 ) (60) 



2J (2vr) 3 ' 



with H(q) being a loop made of the u a propagators: 



[L(p + q/2) + fx 3 ][L(p - q/2) + fi 3 ] [L(p) + ^][L(p) + /i 3 ]J ' v ' 

Eq. (^) defines a peculiar 0(n) sigma model in three dimensions at n = 4. Were n a 
large parameter the fluctuations of the A' field would be suppressed as 1/n. In that case one 
could neglect, to the first approximation in 1/n, the second and higher terms in eq. (|60"D. 
Though n = 4 cannot be considered as particularly 'large' we shall nevertheless leave out the 
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quadratic and higher-order terms in A'. In principle, these terms can be taken into account 
systematically in the 1/n expansion similarly to the way it is done in two-dimensional sigma 



models JTO, M, O] . 



The dependence on the constant part of the Lagrange multiplier fi 3 is given by (see 
eqs.(|6|,|g)) 




n 



(=4) f d 3 p 



log L(p) + /i s 



2 J (2tt) 3 } 
This expression has an extremum in fi determined by the 'gap' equation 



- 2 



d 3 p 



0. 



(62) 



(63) 



f J (27r) 3 L(p) + fj 3 
the integral here being logarithmically divergent since L(p) ^ c|p| 3 at large \p\ (see eq. (|57|)). 



Introducing an ultraviolet cutoff M we get from (p3[) the vacuum expectation value of the 
Lagrange multiplier A = fi 3 where 



fi = const M exp 



57T 



(64) 



9 2 2 (4/c) 

Comparing it with the asymptotic-freedom law (|9]) we see that in order to get the correct 
one-loop Gell-Mann-Low coefficient (at iV c = 2) we need to put 



11 
~3~ 



or 



1 + 



Hp) 



\p\ 
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IT' 



(65) 



We notice that the first Gell-Mann-Low coefficient has been recently calculated in ref. 
H in the Kogan-Kovner ansatz using a different method, and it has been found to be 12/3 
instead of 11/3. This is exactly what we get from eqs.(|>|, |65]) if we deliberately put (as it has 
been done in ||) b(p) = 0. However we emphasize again that there are no reasons to bind 
oneself with such restriction. The actual value of the trial function b(p) at large \p\ should 
come from minimizing the second-order 0(M A g 2 ) term in eq. fllCf ). We have systematically 
neglected such terms above, however they can be systematically collected as well. Those 
terms contain, in particular, the non-abelian commutator terms for the magnetic energy, 
involving 'transverse' gluons. Therefore, it can well happen that the 'best' b(p) —>■ \p\/ll, as 
it should be for getting the correct Gell-Mann-Low coefficient. 

Eq. ( |64"D demonstrates how transmutation of dimensions occurs in the Schrodinger ap- 
proach. We see that it is related to the non-zero vacuum expectation value of the Lagrange 
multiplier field needed to get rid of the unitarity restriction for gauge transformations. This 
phenomenon is similar to what happens in the 0(n) sigma models in two dimensions. Simi- 
larly to those models one needs, strictly speaking, the large n parameter to justify the use of 
the saddle-point method. However, one can hope that n = 4 is sufficiently 'large' so that the 
conclusion that the Lagrange multiplier field A does get a non-zero v.e.v. is not altered by 
fluctuations of the A' field about the saddle-point. These fluctuations can be anyhow treated 
systematically in a formal 1/n expansion. 

It is instructive to separate the perturbative and nonperturbative contributions to the 
vacuum energy by adding and subtracting the perturbative contribution (|48|) . We have 
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2{Nl 



£ gpert gnonpert 
— > 1 

V - V V 
d 3 p \p\ r d 3 p 



(2vr) 3 2 J (2vr 



^3 



(iy2 _ lN L(p) L 2 (p)5logA^(27r) J 



(66) 



y 2p 2 2p 2 SL(p) V 
gnonpert j g a f unc ti ona i of the kernel L(x — y) only (see eq. P2[)), and the variational 
principle comes to minimizing this functional. To perform that one needs to evaluate the 
norm Mu Q%5] ) for arbitrary L(x — y) and then mimimize the r.h.s. of eq. ( pB| ) in L. After the 
transmutation of dimensions takes place, the bare coupling g 2 disappears from all expressions, 
being replaced by a mass parameter u ~ A, see eq. (j64|) . In momentum space L(p, [/,) 
should asymptotically go to (12/11) |p| 3 at \p\ 3> «, however this behaviour cannot follow 
from the minimization of the nonperturbative vacuum energy but should be rather imposed 
'by hands' as presumably following from the minimization of the second-order perturbative 
energy. Apart from this boundary condition, the shape of L(p,p,) is to be extracted from 
the variation of S nonpert which has to be convergent at large momenta. The whole procedure 
seems to be feasible in the framework of the 1/n expansion. 

In the leading order of the 1/n expansion one freezes the Lagrange multiplier A at its 
saddle-point value A = u 3 and neglects the fluctuations of the A' field. The propagator of 
the u a field becomes 

g 2 1 

(u a (p) Uf)(-p)) = —8 aP + ^ 3 . (67) 

Concluding this section we would like to make the following remark. In refs. U 
another method of circumventing the unitarity condition has been suggested. In these refs. 
integration over unitary matrices U has been replaced by 



JDU = J DU^DU* p a 5{\U ] Uf p -5fj = J Da a p j DU^DU* p a expj- J 'dxTi a{U ] U - 1) j . 

(68) 

Here one has N 2 complex fields U, W with 2N 2 real degreees of freedom (neglecting the 
difference between the SU(N C ) and U(N C ) groups), with N 2 restrictions taken care of by 
integrating over N 2 Lagrange multiplier fields a = a'. Taking iV c — ► oo does not help to 
justify any saddle point in a, as it is assumed in refs. §]. Indeed, the effective action 
in a resulting from integrating over the unconstrained fields U, W depends not only on 
N c eigenvalues of varying slowly but also on N 2 — N c — 0(N 2 ) 'angle' variables of 
which fluctuate violently. Meanwhile, the fluctuations of the 'angle' part of influence the 
effective action for the eigenvalues of the a field, and taking N c — > oo does not quench the 
'angle' variables as their number is 0(N 2 ). In short, the method does not work. 

In order to get use of the large N c parameter one has to parametrize the SU(N C ) group 
by introducing 0(N 2 ) 'flat' variables subject to 0(N C ) constraints. In the case of N c = 2 it 
has been explicitly performed above. For larger N c we are, unfortunately, unaware of such 
a parametrization. 

Finally, it should be mentioned that the second derivative of the effective potential ( |6^D 
in A is positive-definite, which may seem to invalidate the use of the saddle-point method. 
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In fact, it is exactly the needed sign of the curvature, since integration over the Lagrange 
multiplier field A in eq. fl56|) actually goes along the imaginary axis. The integration path 
climbs to the maximal value of the integrand (which happens to be exactly on the real axis!) 
and then descends, as it should be in the saddle-point integration. 



9 A queer way to get the linear potential 

Let us now consider two static colour charges in the fundamental representation on top of 
the vacuum state. We shall use the trial wave functional suggested by Zarembo ||. Apart 
from being a functional of the field, it is a function of where we put the probe quark (z\) 
and antiquark (z 2 ); it is also a colour matrix: 



%[A; Zl ,z 2 ] = J DU U"(zi)Up(z 2 ) J DA exp (-^p-) . (69) 



Q\A 
9 

One can immediatelly check that under gauge transformations the wave functional ( 6"9" ) 
satisfies the needed condition (|]). As correctly stressed in ref. || the static potential between 
probe quarks is anyhow a property of the vacuum state and it does not matter much what 
precisely trial functional do we use to 'measure' it, provided it satisfies the Gauss condition 
(H). Eq. (|69"D is a simple and natural choice. 

By definition, the potential between static sources is the hamiltonian averaged with the 
wave functional (^)9]) minus the vacuum energy: 

v( _ ^(z u z 2 )\H\^(z 1 ,z 2 )) (* \H\9 ) 

1 1 2) (y(z l7 z 2 )\y(z u z 2 )} (tf |*o> • 1 ] 

The norm of the ^l{z\^ z 2 ) state is H 



M( Zl -z 2 ) = (V(z u z 2 )\^(z u z 2 )) 



DU / DA Tr 



U( Zl )U\z 2 ) 



exp 



Q[A\ Q[A L 



9" 



9" 



(71) 

where we have again, as in the vacuum case, introduced the relative gauge transformation 
U = U\U 2 . Performing the same one-loop approximations as in sections 5 and 6, it is easy 



to see that the contributions from the Ai field get cancelled in eq. (|70|) , and we are left with 
the correlation functions of the sigma model. Let us introduce a short-hand notation for 
these correlation functions: 



«<?» 



JDU exp (-^JJdx dy L(x - y) Tr W(x)U(y)) O 



JDUexp (-^JfdxdyL(x - y) Tr W(x)U(y 
The static potential can be written as 



(72) 



V{ Zl -z 2 ) = - 



JJdx dy Sjx - y) ((Tr U{ Zl )W(z 2 )Tr W(x)U(y))) 
<r £ L {(TrU( Zl )W(z 2 ))) 

ffdxdyS(x-y)((TrW(x)U(y))) 



(73) 
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In the leading order of the 1/n expansion the integrals in eqs.([72|, |73|) become gaussian with 
the propagator given by eq. (^)- There are three possible contractions of the u a fields in 
the numerator of eq. (|73|): one of them gets cancelled by vacuum subtraction in (|7^); the 
other two are equivalent. Using the u a propagator (|67D we obtain 



E{r) AJV _ N o f d 3 p i[p . r) 1 _ 2 [°° j„ sin(pr) 



^( r ) = T77T, N(r) = 2 — ^e %{ P- r > - - - = — / dp 



N(r) } KJ 7 (2vr) 3 L(p) + /i 3 2 Wo F p L(p) + fi 



,3 ' 



Let us first investigate the short- distance behaviour of V(r). Small r correspond to large 
p where L(p) = (12/ll)|p| 3 according to eq. (|65|). We have at small r 

E(r) -> — , Mr) -> — 4 In -, V(r) -> — Jjir. (75) 

v ; 4vrr' v 7 12 7T 2 r' 17 4vrr f In i V ; 

3 r 

This is almost the correct Coulomb potential modified by the running of the gauge coupling 
constant. The only difference with the true asymptotics at small r is the absence of the 
(N 2 — l)/(2iV c ) = 3/4 colour factor in eq. flT5|). The difference 1 — 3/4 = 1/4 is due to our 



use of the leading order of the 1/n expansion at n = 4; simultaneously it gives the idea of 
the numerical error. 

The behaviour of V(r) at large separations depends crucially on the analytic properties 
of L(p). The asymptotics of the Fourier transforms (|74|) are defined by the singularities of 
the integrands in p 2 . If L(p, //) + /i 3 has a singularity at p 2 = (for example, is odd in |p|) 
the large-r asymptotics of both E(r) and N(r) are determined by the integration origin at 
\p\ = 0: both quantities exhibit a power decay at large r. Examples we have considered all 
show that their ratio V(r) is also decreasing as a power of r. This is not what we want for 
confinement. 

If L(p, //) + /i 3 is even in \p\ one can expand the integration range in \p\ in eq. (|74|) 
from (0, oo) to (— oo, oo). The asymptotics will be then determined by the positions of the 
singularities of the integrands in the \p\ plane. If these singularities have a real part both 
E{r) and N(r) and consequently V(r) are oscillating. Probably, this is also unacceptable. 

If L(p, fi) + /i 3 has zeros on the imaginary axis of \p\, say at \p\ = ±m, E(r) and N(r) 
decay exponentially as ~ exp(— rer). The key point is that if the integrand of N(r) has 
a single pole at \p\ = ±i/c, the integrand of E(r) has a double pole. Taking the residue 
of the double pole one can differentiate exp(z|p|r) and get an extra factor of r in E(r) as 
compared to N(r). The same mathematical mechanism works if L(p, p) + fi 3 has a cut at 
some imaginary value of \p\. This is the way one obtains an infinitely rising linear potential. 

To give an example we choose 



T , n. 3 12 
L[p,fi)+fi = — 



2 l U V 2 



P +1^) A* 



(76) 



satisfying the above considerations and also the asymptotics (12/ll)|p| 3 . The coefficient 
inside the square brackets has been chosen so that b(p) is not singular at p = 0. The 
corresponding potential V(r) is plotted in Fig.l. 
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Figure 1: Static potential following from eq. (|7B]). The axes are in units of fi obtained from 
the transmutation of dimensions. The standard string tension a ~ (440 MeV) 2 corresponds 
to n « 600 MeV. 

Both the Coulomb part fl75|) and the linear part are clearly seen. In mathematical terms, 
both Eir) and N(r) decay exponentially as ~ exp(— const fi r) but N(r) is additionally 
suppressed as 1/r. Being inversed, it gives a linear rising potential. 

Let us put it in a more general way: In the original formula for the potential (|70|) both 
the numerator and the denominator have to decay exponentially if there is a mass gap in 
the theory. However the denominator, that is the norm of the state with colour sources 
(^(^i, z 2 )\^(zi, z 2 )), has to fall off as a factor 1/r faster than the average of the hamiltonian, 
($(zi, z 2 )\'H\ty(zi, z 2 )). Of course, one can always artificially normalize the wave functional 
to unity but that does not seem to be a natural thing to do. 

We have called this way of getting confinement 'queer' because it is very different from the 
usual Euclidean standpoint where one has to demonstrate a tremendously small exp (—Area) 
behaviour of the Wilson loop or a similar behaviour of the Polyakov's loops correlator. This 
is very difficult both theoretically and in practice. However in the Schrodinger approach 
getting the linear potential is not at all queer but rather quite natural. 



10 Adjoint sources 

Let us first of all introduce the gauge transformation in the adjoint representation: 



O a \U) = 2Tr \uH a Ut l 



QacQbc = QcaQ 



ca/^icb ^ab 



(77) 



We take the wave functional for a state with two adjoint sources sitting at points z\$, in 
the form similar to eq. (|69|): 



V ab [A; Zl ,z 2 ] = J DU O ac { Zl )O bc {z 2 ) jDAexp 
It satisfies the needed gauge transformation law, 



Q[A 



(78) 
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^ ab [A s ; Zl ,z 2 ] = 0-\S( Zl )) - ^' b '[A- z u z 2 ] [0{S{z 2 ))\ 



,a'b' 



ib'b 



(79) 



Proceeding in the same fashion as in the previous section we arrive to the following 
expression for the potential for static adjoint charges: 



V adj { Zl -z 2 ) = -- 



JJdx dy S(x - y) (( hi U^ Zl )U(z 2 )Tr U{ Zl )U\z 2 ) - ll Tr U^(x)U(y))) 



((TtW( Zi )U( Z2 )TtU( Zi )W( Z2 ) - 1)) 



ffdxdyS(x-y)((TrW(x)U(y))) 

((D) 



10) 



Passing from integration over unitary matrices U to the flat variables u a and the Lagrange 
multiplier A one observes that there are much more possible contractions of u a fields than 
in the case of fundamental sources. Not going into details here we mention that thanks 
to the 'gap' equation contributions potentially possessing the confining behaviour get 
cancelled. There is no linear potential between static adjoint sources even if there is such 
for the fundamental charges. 



11 What is to be done next 

We believe that the variational approach may become a powerful alternative to the Euclidean 
study of gluodynamics. To put it on a more solid basis one needs: 

• To compute within the variational ansatz given by eqs.(|l^, O, [14]) the first perturbative 
correction to the vacuum energy, which is of the order of g 2 M 4 . Only the asymptotics of 
the trial functions a(p), b(p) at large \p\ enter this calculation. The function a(p) = \p\ 
is determined from the first-order calculation (see above), the function b(p) is not. 
Computing the 0(g 2 ) correction will fix this function. There is a good chance that the 
best b(p) will turn out to be |p|/ll as required by asymptotic freedom, eq. (|65|). 



• To compute 1/n corrections to the saddle-point approximation to the sigma model (p6|) . 
It will enable one to learn, by minimizing the r.h.s. of eq. (|66|), the nonperturbative 
propagator L(p, p) + fi 3 . 

Both computations seem to be feasible. There are also obvious applications of the vari- 
ational approach: 

• The glueball masses and wave functions 

• The 6 dependence of the vacuum energy, or the topological susceptibility of the vacuum 

• The area behaviour of the spatial Wilson loop (?) 
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12 Conclusions 



The gaussian variational ansatz for the vacuum wave functional is, on one hand, sufficiently 
simple to allow its study by conventional field-theoretic methods but on the other hand, it 
is probably sufficiently 'rich' to incorporate many desired features of the Yang-Mills theory. 

We have shown that by including the longitudinal structure in the gaussian kernel one 
can, in principle, get the correct one-loop /3-function of the Yang-Mills theory within the 
variational approach. 

The famous transmutation of dimensions arises very natural in the Schrodinger picture. 
It results from a phenomenon similar to what is well known in 0(n) sigma models in two 
dimensions, namely from the appearance of a nonzero vacuum expectation value of the 
Lagrange multiplier field needed to get rid of the unitarity constraint on the gauge transfor- 
mations. We have performed it explicitly for the 577(2) gauge group; unfortunately we are 
unable to generalize the method to higher N c . 

The appearance of the mass gap through the transmutation of dimensions is good news, 
however it is not sufficient by itself to get a linear potential between source quarks. The 
confinement requirement can be formulated as a simple analyticity property of the propagator 
in the sigma model, once the Lagrange multiplier gets a nonzero v.e.v. This propagator can 
be systematically studied within the 1/n expansion, again similar to what has been developed 
in d = 2 models. 

The variational approach allows certain freedom in action. Instead of working hard trying 
to establish the best parameters / functions within a variational ansatz, one can choose them 
to be what one likes, for example, taking the propagator in the form of eq. fl76|) leading to 
a beautiful static potential shown in Fig.l. After all, the variational principle says that 
one gets an upper bound for the energy for any trial function. Using, for example, eq. (|76|) 
it is interesting to check if it leads to reasonable glueball states, reasonable topological 
susceptibility, and so on. 

The most challenging theoretical question is whether the variational wave functional gives 
rise to the area law for the spatial Wilson loop. To check it, the representation of the Wilson 
loop via the nonabelian Stokes theorem of ref. ]D|] would be probably helpful. 
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